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In the flurry of experiments looking for topological insulator materials, it has been recently discovered that 
some bulk metals very close to topological insulator electronic states, support the same topological surface states 
that are the defining characteristic of the topological insulator. First observed in spin-polarized ARPES in Sb 
(D. Hsieh et al. Science 323, 919 (2009)), the helical surface states in the metallic systems appear to be robust 
' to at least mild disorder. We present here a theoretical investigation of the nature of these "helical metals" - 

, bulk metals with helical surface states. We explore how the surface and bulk states can mix, in both clean and 

disordered systems. Using the Fano model, we discover that in a clean system, the helical surface states are 
not simply absorbed by hybridization with a non-topological parasitic metallic band. Instead, they are pushed 
away from overlapping in momentum and energy with the bulk states, leaving behind a finite-lifetime surface 
resonance in the bulk energy band. Furthermore, the hybridization may lead in some cases to multiplied surface 
state bands, in all cases retaining the helical characteristic. Weak disorder leads to very similar effects - surface 
states are pushed away from the energy bandwidth of the bulk, leaving behind a finite-lifetime surface resonance 
in place of the original surface states. 
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I. INTRODUCTION 



, The prediction and subsequent discovery of bulk topological insulators (TI) has galvanized the condensed matter community, 
not in the least because of the unique physics of the protected edge states on their surfacai— . Many materials adjacent to 
TPs, however, are bulk metals with the same helical surface states, e.g., Sb^^-^^, Bio.giSbg.og^^^, Bi2-2;Mn2:Te3-'^, and even 
undoped Bi2Se3-'-''^'*, which initially was thought to be a wide band-gap TI, but now seems to have a small but measurable bulk 
^ Fermi surface. In such materials one would naively expect that a bulk Fermi surface would simply swallow the surface states. 
O Nevertheless, the same helical surface states of the TI phase, continue to appear in them, now coexisting with a bulk Fermi 
I ^ I surface. We suspect that more examples will emerge, especially since the helical surface states are less sensitive to disorder than 
ordinary (non-helical) surface states (of the likes of those found on the (111) surfaces of Cu), as has been already suggested by 
CO surface-doping Sb with Pottasium^**. 

The prevalence of metallic systems with surviving helical surface states led us to ask: what exactly happens when a surface- 
state, which is the result of a TI bulk, is allowed to hybridize with a 'parasitic' non-topological metallic band. As we shall see 
below, not only do the helical edge states survive, but they can actually multiply. When mixed, the surface states are simply 
. pushed away from overlapping in energy and momentum with the bulk states. If the surface and bulk do not overlap to begin with, 
the mixing is not effective, and the surface states are only slightly modified. In those areas of energy and momentum overlap 
between the bulk and surface states, new "exiled" surface states appear above and below the confines of the metallic band, and 
in place of the original surface state a "ghost" surface resonance remains, with a finite (and often very short) lifetime. Thus in 
one momentum value, we may see ARPES signatures of two surface states at energies above, and below, the metallic energy 
!• range (see Fig.[Tll, in addition to a surface resonance within the bulk metallic band. The exiled states, as well as the remaining 
finite-life-time ghost resonances will retain the odd-number distinction of surface bands, characteristic of helical surface states, 
comprising an odd number of Kramer's pairs of bands in 2D, and an odd number of Dirac cones in 3D. 
5—1 ' In this manuscript, we explore the metal vs. edge-state struggle by first constructing 2d examples where a metallic band 
appears at the same energy as an edge state, and numerically investigating its emerging spectral structure. Next, we approach 
the problem analytically by constructing a generic model for a helical surface state interacting with a bulk metal, based on the 
Fano model^. From both approaches, the generic picture of a helical metal arises, as is summarized in Fig. [T] 



n. CONSTRUCTION OF 2D HELICAL METALS 



We begin our study by attempting to construct disorder free theoretical models in 2d where helical edge states coexist with 
a bulk metal in both momentum and energy. Such helical metals can be achieved in at least two different ways. One possible 
construction is to take a model of a TI, add a new (initially decoupled) partially filled band, an then mix it with the TI bands. A 
second possible construction is to add a momentum-dependent chemical potential without adding additional degrees of freedom, 
such that the gap closes somewhere in the Brillouin zone (BZ). This can always be achieved in a tight binding model by adding 
appropriate hopping terms that produce purely diagonal terms in the multiband hopping model. We will consider an explicit 
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FIG. 1 : An energy-momentum sketch of the fate of surface states when they overlap with a parasitic metallic bulk band. The solid (light green) 
regions mark the bulk density of states, the top and bottom bands being bulk states of a topological insulator, and a parasitic metallic band 
in the middle. The original surface state branch (dashed line) overlaps with the middle (metallic) band. The surface state branch hybridizes 
with the bulk, leaving a diffuse resonance in the bulk band energy range, as well as two sharp edge-state branches appearing above and below 
the parasitic metal energies. This demonstrates that the topological surface states, rather than being obliterated, are simply pushed away from 
overlapping with the bulk. 




FIG. 2: A two dimensional toy model for the investigation of a helical metal. The model is based on the Kane-Mele model'"'' on the honeycomb 
lattice, which is denoted by dashed lines. The Kane-Mele model consists of nearest-neighbor hopping ti (green) between the two (a, h) 
sublattices of the honeycomb lattice (denoted by empty and filled circles respectively), and a complex second neighbor hopping (a spin orbit 
coupling term) with opposite sign when clockwise (— iA, blue) and counterclockwise (+iA, red), as indicated by the curved arrows in the figure. 
In addition to sites of the honeycomb lattice, we include a new set of sites (c) at the centers of the hexagonal plaquettes of the honeycomb 
lattice, denoted by cross marks. The c-sites form a triangular sublattice, and nearest neighbor hopping between them t2 (brown) forms a 
metallic band. To explore the interplay between the topological insulators helical surface states and the bulk metallic band we mix the two 
systems by allowing hopping between the c-sites and the honeycomb lattice sites, is (purple). 



example of the former, demonstrating the construction in 2d. An explicit example of the latter construction is relegated to the 
supplementary material (see Section IaTI i. as the main physical features of it are no different than in the one example we show 
here). Our 2d examples can easily be generalized for models of topological insulators in any dimension. 

Our starting point is the first 2D model of a topological insulator, the Kane-Mele model'-^, defined on a honeycomb lattice. 
The first construction proposed above can be realized by considering this model coupled to a very simple metallic band: we add 
lattice sites at the centers of the honeycomb plaquettes, as shown in Fig. |2l so that the new sites are now those of a triangular 
lattice. To form a metallic band we allow hopping between the nearest-neighboring new (triangular) sites. With a Fermi energy 
crossing anywhere in this band, we have a single Fermi surface, centered about q = 0. We then couple between the two parts 
of our model, by allowing hopping between the honeycomb lattice sites and the new triangular lattice sites (see Fig. |2]l. The 
Hamiltonian is 



^3 E 



(ij)a 



h.c. 



Aa^jP^aP^ij + '^■C- + (a — > 6) 



(1) 



where i,j denote all the combined lattice sites, a, &, c denote the fermion operators on the three sublattices (a, h for the honey- 
comb, and c for the triangular metal, see Fig.|2|i. In addition, spin indices are denoted by a, (3, and Vij is as defined in Ref llO 
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(e)Zigzag edge, (3 = 0.3 (f)Armchair edge, = 0.3 



FIG. 3: Spectrum of the modified Kane-Mele model in Eq. (TJ, with a metallic (triangular) band added (Fig. Here we use the parameters 
ti = 1, X = 0.5, t2 = 0.2, h = 0.4, and vary ts. We show spectra for both zigzag l |3(a)l[3(c)l[3(e)l l and armchair edge strips l |3(b)l[3(d)][3fflt , 
with mixing parameter values is — 0, 0.03, 0.3. The strips are finite in one direction, and periodic in the other, so the lattice is wrapped around 
a cylinder. The edge-state branches on different sides of the sample are denoted by thick lines, dark (blue) and light (orange) respectively. All 
bulk states are colored light gray. Helical edge state branches can still appear (at least when the Fermi energy is shifted), as can be easily seen 
by counting the (odd) number of surface state branches on one side of the strip crossing E = between q = and the BZ midpoint (q = n 
for zigzag, q = 7r/\/3 for armchair). 

(and shown in Fig.|2]i. For convenience we have included an independent chemical potential h for the triangular lattice. 

Following closely the procedure in Refs.ljllsl we calculate the spectrum of Eq. ^ using exact diagonalization of strips that are 
terminated either at a zigzag or an armchair edge, as these are representative edge cuts of generic boundaries of the honeycomb 
lattice. We calculate for one spin polarization in order to avoid clutter, as is a good quantum number in our model, and 
deviations from this are immaterial to the physics we explore here, and identify explicitly the eigenstates localized at each edge 
of the finite strip. We present the spectrum of the model in Fig. |3] when the honeycomb-triangular hybridization parameter is 
zero, weak and strong (t^ — 0, 0.03, 0.3 respectively), while keeping the other model parameters fixed at ii = 1, A = 0.5, t2 = 
0.2, h = 0.4. We choose the Bravais lattice vectors to be ai = ax, a2 = | [—x + y/3y). The strip widths we use are 30 unit 
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cells, and are terminated in a symmetric fashion. For the zigzag geometry, states identified as edge states on one side of the 
system have at least 0.85 of their total weight within a distance of y = 30| from the lower edge (at y = 0), and the edge states 

at the other side of the sample have at least 0.85 of their weight fraction above y — 30 "^^^ (for the armchair edge strip, we 
take the limits for this procedure to be a; < 30^-^ and x > 30-^, and consider 0.7 of the total wavefunction weight rather than 
0.85). 

a. Weak hybridization. In Figs. |3(c)[ |3(d)[ we plot the spectrum of bulk and edge state for an armchair and zig-zag strip 
with weak hybridization, — 0.03. Comparing with the decoupled spectrum in Figs. |3(a)]|3(b)[ the numerical results clearly 
show that essentially the same helical surface states appear in our model without a bulk gap present - the helical surface states 
are only weakly perturbed. Only in those areas where the helical surface states overlap in both energy and momentum with bulk 
states, do we see any appreciable change. Since different momenta cannot mix in the absence of disorder, it is not surprising the 
surface states are unaffected if they do not overlap in both energy and momentum. 

b. Strong hybridization. We consider also strong coupling (t^ = 0.3) between the honeycomb and triangular subsystems. 
The spectra, shown in Figs. |3(e)|3(f)| in which, surprisingly, the surface state branches have been pushed out of those regions 
where the bulk states reside, away from their original position, toward those regions of the momentum-energy diagram where 
bulk states are absent. This results in the doubling of the number of surface state branches for a single momentum, as can be 
seen in Fig. |3(e)| in the range 2.4 < q < 3.2. One branch spans the region below the metallic density of states, and another 
above. 

It is noteworthy that for other parameter choices in this model, one can also find instances where the number of co-moving 
surface states changes upon bulk-edge coupling, but still with an odd number of surface branches on each side of the strip (in the 
supplementary material, Section lA2l an example is given where the number of surface states changes from 1 to 3 - still retaining 
the odd number of branches). 

Our numerical results suggest a tendency of surface states to be pushed away from overlapping in (surface-projected) momen- 
tum and energy with parasitic bulk states originating from a non-topological band. Not only do they seem to persist, in some 
momentum number they seem to multiply. Motivated by this observation, we turn next to an analytic treatment of the mixing 
between surface states and bulk states. 



III. THE FANO MODEL APPLIED TO EDGE-STATES - BULK MIXING 



Our numerical results confirm that edge states have remarkable resilience. Even in the absence of a band gap in the material, 
precise edge states form in empty spaces in the projected energy-momentum state diagrams. Let us now approach the problem 
in its idealized form analytically. For this purpose we write a model consisting of a surface state branch and a decoupled bulk 
metallic band, and then allow them to mix. This most certainly describes the specific model we have discussed here, since an 
effective low-energy continuum theory would encompass precisely these elements.'*^. We label the transverse momentum by q, 
and the momentum in the direction perpendicular to the wall by k. In the clean limit, the bulk and surface states are mixed by 
hopping matrix elements which preserve q, and therefore states with different q do not mix. This allows us to treat individual 
edge states separately. 

The model outlined above takes the form of the well known Fano model^. Making use of the path-integral formulation, the 
action for the model is 



[iw„ - e(q)] + x^xl^^n - E{k,q)] + / [i:''xgik,q) + h.c 



(2) 



where the fermionic Matsubara frequencies are w„ = tt ^"^''" , and the surface and bulk states are denoted by the Grassmann 
fields tp = ^), and x — xi'i^n, k, q) respectively, their dependence on momentum and Matsubara frequency suppressed 

for the sake of brevity. The surface and bulk state energies are e(q) and E{k,q) respectively, and the coupling between them 
is g{k, q). The exact solution of the Fano modeP*", can be most easily achieved by integrating out the bulk degrees of freedom. 
Since the action is quadratic, this can be done exactly, resulting in an effective action for the surface state degrees of freedom 
(-0) alone: 



l5(^,q) 



i?(fc,q) 



(3) 



This action reveals the fate of the spectrum of the surface states. The retarded Green's function extracted from the effective 
action above is 



^ret 



Jk 



l.9(fc,q)l= 



E{k,q)+iS 



(4) 
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FIG. 4: Typical spectral weight function of the Fano model. We take the spectral weight in Eq (O, with parameters Ei = —1, E2 = +1, 
the surface state energy e = 0.5 (within the bandwidth). For uj values outside the bulk bandwidth we take A = 10~* to approximate the 
delta function outside the bandwidth. For the continuous (red) curve we take \g\'^iy{EF) = 0.04, and for the dashed (blue) curve we take 
Igl^f^Ep) = 0.2. At the lower coupling value (continuous curve), the widened peak near uj — 0.5 is the remnant of the original surface 
state, and the modest weight in one delta function already appears very near the band edge (the weight of the other delta function is too small 
to discern). For the stronger coupling value (dashed line) traces of the surface state are largely obscured in the bulk bandwidth, but both the 
delta-functions outside the bandwidth appear prominently. 



From which the spectral function ^(q, w) = —2Im [Gret{^.^ ^)] can be extracted. 

The most important feature of the spectral function A(q, uj) is that the original surface states, when their energy is within the 
bandwidth of the bulk, acquire a lifetime, which is roughly r ~ \g\iu{Ep) (where ^{E) is the density of states, and Ep the Fermi 
energy). This lifetime describes the typical time scale in which the edge probability density leaks into the bulk states. In addition, 
new delta functions appear at energies outside the bandwidth of the bulk metallic band both above and below. We illustrate this 
by considering a constant coupling g, and a uniform density of states v{Ep), in a bulk band with energies Ei < E < E2- The 
spectral function is 



(5) 



where A — \g\'^i'{Ep), and /to(S) = — ttA for Ei < lu < E2, and zero if lu lies outside the bandwidth of the bulk band, thus 
producing the aforementioned delta functions. Fig. |4] shows the spectral function for particular values of the parameters. If the 
original surface state energy lies outside the bulk bandwidth, it remains as a delta function outside the bandwidth, slightly shifted 
from the original surface state energy, and with most of its spectral weight retained. In addition, a new delta-function peak will 
appear on the other side (in energy) of the bulk band, with a small spectral weight, and some spectral weight will appear within 
the bulk bandwidth. 

These features of the Fano model depend only weakly on the details of the density of states, as can be demonstrated by 
considering a more realistic bulk band dispersion E{k) — — /i — 2t cos(fc). Using the dimensionless variable u = {uj + n)/2t, 

the integral for the self-energy can be solved analytically, yielding i?e(S) — V^S^ ^'^'^ ImiYi) = when \u\ > 1, and 

2 

Re{T,) — and Im(Jl) — — ^-^£== when \u\ < 1. The resulting spectral function is qualitatively no different than that of 
Fig.m which corresponds to Eq. Q. 

The remarkable structure the Fano model implies can be further intuitively understood by considering an impurity state (d) 
coupled to a flat band of bulk states (/„) at zero energy, = ed^d + .gX^n [fn'^ + h.c] (the coupling g is taken constant 
without loss of generality), with n — 1 . . .N. The characteristic polynomial of the Hamiltonian matrix Det{E — H) = 
{—1)^^^E^^^ (£'^ — eE — Ng"^), has — 1 zero modes left from the original N flat band states, and two roots at i? = 

e/2 ± outside the flat band, which, assuming Ng^ ^ e^, become E w ±^/N g + e/2. As in the Fano model 

above, exact energy eigenstates appear above and below the band. While this is an extremely artificial example (flat band states 
being most sensitive to coupling to other states), it demonstrates how the salient features of the Fano model can be understood 
- the impurity state "hijacks" one effective mode from the band, and mixes with it, producing two eigenstates with energies 
outside the band, which correspond to reduced-weight delta functions in the spectral function for d. 
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A. Fano model implications 



We can interpret the numerical results of Section |II]in terms of the salient features of the Fano model as follows. The Fano 
model tells us that for those edge-parallel momenta q where the surface states and bulk states do not overlap in energy, the 
surface states are only slightly shifted in energy. This is demonstrated in those regions in Fig. [3] where the original surface states 
did not overlap with the bulk states. On the other hand, if the surface states overlap in energy with the bulk state energies, the 
surface states are "exiled" from the bulk bandwidth, and form states at (very different) lower and higher energies, in addition to 
leaving "ghost" surface resonances overlapping with the bulk states, close to where they originated. The "exiled" states are still 
surface states, mixed with a superposition of bulk states localized at the surface, since any eigenstate existing in a bulk gap in 
energy-momentum space, must be an evanescent wave into the bulk. The surface state spectrum e(q) is a continuous function 
of q, and assuming the coupling g{k, q) is also a continuous function of momentum, the exiled surface states will also form a 
q - continuous energy branch. The helical nature of the surface states will also be preserved in the new exiled surface states: it 
will correspond to an odd number of Dirac cones in 3D (and an odd number of Kramer's pairs of bands in 2D). The exiled states 
are evident in the numerics in Fig. |3(e)| and Fig. |3(f)| in particular in Fig. |3(e)| in the range 2.4 < q < 3.2. Indeed, the exiled 
states retain the surface state branch continuity, as well as their helical nature. The ghost surface resonances are not pure surface 
states, and so are not distinguished in the spectra in Fig.|3] In order to identify the "ghost state" signatures in our numerical, we 
probe the explicit spectral weight function of our model in the next subsection. 

Experimentally, we expect that on different facets of the crystal the two scenarios could be realized, and so if a helical surface 
state appears on one facet of the metal, crossing the Fermi energy, but not on the other facets, it may simply be significantly 
shifted in energy, and can in principle still be observed outside the energy range of the bulk band. 



B. Numerical evidence for the surface resonance 



In this section we will present numerical evidence for the "ghost state" surface resonance in the model described in section Ull 
First, we calculate explicitly the surface state spectral weight function Q described above in the Fano model analysis. This is 
done by calculating a probability distribution to find the original surface state eigenstates at a given energy and momentum in the 
new, coupled spectrum. We take the surface state eigenstates V'o(q) in the decoupled case, and calculate their overlap with the 
various eigenstates of the coupled system, squared. In particular, if we denote by 0^ (q) all the eigenstates of the coupled system, 
with energies uj, the spectral weight is A(q, = | (V'o (q) 1 4''^ (q)) I ^ ■ The results of our numerics, for the surface state branch 
along one edge of the sample in the strong coupling case (^3 = 0.3), are shown in Figs. |5] A faint but discernible diffuse peak 
is seen to overlap with the bulk states in between the upper and lower exiled surface state lines, thus confirming the predictions 
from the Fano model analysis. 

The spectral weight function is not necessarily what ARPES or STM will measure. In order to give a clear experimental 
signature that can be measured, we perform one additional numerical analysis. We take the density profile of each eigenstate 
(/)(^(q), and convolve it with a weight factor f{y) = e^^^, to yield the total weight J^{io, q) = |(0cj(q)|/(2/)|0Lj(q))P- The 
decaying exponential mimics the finite penetration depth surface probes can achieve. Advantageously, the calculation of J', is 
unbiased by the "band archeology" in calculating yl(q, w), which required comparison with the eigenstates of the decoupled 
system. We plot bfq) versus ui and q in Fig.|6l Those eigenstates with a significant part localized at the surface should 
have a sizable value of J'. Indeed, we see in Fig. |6] that between the exiled surface states the signature of the ghost surface 
resonance appears. In conclusion, careful analysis of our numerical results shows that surface state probes could identify the 
"ghost" surface resonance. 



IV. DISORDERED HELICAL METAL WITHIN THE FANO MODEL 



Going beyond the clean limit, we now consider edge-bulk hybridization in the presence of disorder A random disorder 
potential scatters surface states into other surface states, as well as into the bulk. For simplicity we ignore, however, disorder 
scattering between bulk states, concentrating only on the fate of the surface states. Within the weak-disorder approximation^ 
we will see that the Fano picture by and large still applies. 

Disorder scattering can easily be quantified by slightly modifying Eq. ©I 



[lu^ - e(q)] + / x^x - E{k, q)] + / [xH' 9{k, q, q')* + h.c] + / V^^'^(q, q') 



(6) 



where we use the shorthand notation il^' = ip{iujn, q') (and similarly for ip" later on). The coupling term g{k, q, q') no longer 
conserves transverse momentum, and is determined by the overlap of the surface and bulk (x) state wavefunctions repre- 
sented by the various fields weighted by the random potential, g{k, q, q') = V'q' (r)* J/(i")Xfe,q(i")- Integrating out the bulk 
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FIG. 5: Numerical calculation of the surface state spectral function. We take the surface state eigenstates in the decoupled case (i/)o(q)) and 
find their distribution among the eigenstates of the coupled system (?!)tj(q) by calculating the overlaps A{q,Lo) — |{i/'o(q)|</>tj(q))P- This 
is expected to approximate Eq. l|5j. We use the same parameters as the strongly hybridized zigzag edge strip in Fig. |3(e)| and the overall 
spectrum is manifestly identical in the two images. We plot the surface state spectral weight function versus transverse momentum (horizontal 
axis) and energy (vertical axis). Dark (blue) points have A{uj) > 0.06, very light (green) points have 0.06 > A{uj) > 0.02, and all points 
with 0.02 > are gray. The exiled surface state branches are manifest, and the ghost spectral resonance is identifiable as a faint peak 

overlapping with the bulk band, between the two exile branches. In addition, the lower spectral weight points map out the shadow of the full 
bulk spectrum. 




FIG. 6: The 'ghost' surface resonance observed in surface probes. The quantity J measures the weight of each eigenstate (jj^j convolved with an 
exponentially decaying weight factor, mimicking the penetration depth of the experimental surface probes J'{uj,q) = j {(f>^ {ci)\f{y)\(f>^{q))\'^. 
Here we plot J versus transverse momentum q and energy uj, both in a 3D (a) and 2D (b) plot. The gray-scale color scheme is such that higher 
values are colored dark gray. We find significant values of J only for the surface states and the surface resonances, which are manifestly 
positioned in between the exiled surface states (most clearly in the 2D plot (b). In addition, the very light gray points map out the bulk states. 
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fermions, we are left with the effective action: 




J2 I N« - e(q)] + / ^li>q'V{ci, q') - / 





(7) 



Next we average over disorder, assuming the disorder has a Gaussian distribution U{y)U{y') — C(r — r') = ii?5{y — r'), 
and surface and bulk wavefunctions ^ g«q R-kl/d ^jj^j ^ ^iq-n+ikx (^vjjgj-g jjjg edge-parallel and perpendicular coordinates 
have been separated as r = {x, R), and d is the skin depth of the surface states). Finally, expanding in the limit of weak disorder 
we find the leading contribution from disorder is a term quadratic in ip (higher order terms in tp are weaker in this limit) 



where u differs from u by some numerical constant. It is important to note that this leading term originates solely from the bulk- 
surface state scattering, and can also be derived by simply taking the disorder mean of the action in O. This gives a Green's 
function nearly identical to that of dU, but we notice that while transverse momentum is now once again a good quantum number, 
disorder, even after averaging, couples all bulk states to any one of the surface states. Therefore, the separation in energy and 
momentum picture no longer holds here, and any surface state that does not exist in a full energy band gap will suffer the effects 
of a Fano mode impurity state in the energy bandwidth of the bulk band. From the original surface state branch, a "ghost" feature 
will remain in the spectral function, and "exiled" surface state delta-functions appear outside the bulk bandwidth. When disorder 
is very weak, the "ghost" surface states will actually be the most noticeable feature, and will simply seem as ordinary surface 
states with a lifetime to leak into the bulk, determined by the strength of the disorder 

A simple experimental test for our predictions follows from the observation that the scattering strength is proportional to 
impurity density: cx pimp- In samples where the disorder scattering of surface states into the bulk dominates over other 
effects (phonon scattering and electron-electron interactions), the surface-state lifetime must be inversely proportional to the 
impurity density. Very recent work^*^ has already measured a quasiparticle lifetime for the surface states in Bi2Se3, which is 
one of the main candidates to be a helical metal as mentioned above, and concluded from the energy dispersion of the surface 
state lifetime that disorder scattering seems to be the dominating scattering mechanism. It would be particularly interesting in 
such experiments to search for the exiled states, which should be above and below the bandwidth of the bulk band. We leave a 
theoretical study of the nature of the exiled states in disordered helical metals to future work. 

The difference between helical surface states and non-helical surface states becomes paramount when considering disorder 
Non-helical surface states, in the clean limit, consist of an even number of bands that can mix via time-reversal preserving terms. 
In 2D materials with ID surface states, non-magnetic disorder, of the likes we consider here, can back-scatter between two 
such bands since their spin configurations are non-orthogonal, and therefore strongly localize them (although a strong spin-orbit 
interaction may mitigate this effect^). A helical edge state on the other hand is a chiral ID conductor, which suffers no backscat- 
tering. In 3D materials the surface states are 2-dimensional, and non-magnetic disorder would lead to weak localization"*". 3D 
helical surface states, however, always exhibit anti-localization of a single Dirac cone. Non-helical surface states may also be 
in an anti-localization class due to spin-orbit coupling, but unlike the helical states, they are (topologically) smoothly connected 
to a spin-rotationally symmetric 2DEG, which suffers weak localization. The additional protection that helical surface states 
exhibit against localization effects indicates that the surface-bulk Fano effects will dominate over localization effects, whereas 
the opposite may be true for non-helical systems. 



Topological band insulators are characterized by their surface state properties. Helical states on the surface of a TI have odd, 
rather than even, number of either Kramer's pairs of Id surface channels in a 2d material, or 2d surface metals with a Dirac 
dispersion in a 3d material'"-'''^ '^. In our work we demonstrated that even materials that are generically metals retain features of 
the helical structure of the surface states, and therefore a similar classification can apply to metallic states. 

Furthermore, we find that the mixing between surface states and bulk states in metals results in a generic rearrangement of the 
surface-state spectrum. At energies where surface and bulk states overlap and mix (either when they overlap also in momentum 
parallel to the surface, or due to disorder) the surface states are reduced into surface resonance "ghosts", with a diffused spectral 
function peak centered near the original energy of the surface states. In addition, the bulk-surface mixing produces "exiled 
surface states" outside the confines of the bulk parasitic metal energy bandwidth. Therefore, the surface-state spectral structure 
emerging from this mixing, contains for each surface momentum q (where surface states exist) two sharp surface states above 
and below the parasitic metal bands where there is a bulk gap (for that momentum number, or, when weak disorder is at play, a 
complete bulk gap), and a diffuse surface resonance overlapping with the metallic density of states. Surface resonances, despite 




(8) 



V. CONCLUSIONS 
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being immersed in the bulk states, can still be identified in ARPES measurements due to the fact that they have very weak fcj^- 
dispersion"*' (and ideally none), in contrast to bulk states. Following the surface resonance peak for different q's should roughly 
parallel the sharp surface states spectrum above and below the band. The location of the sharp exiled states and the diffuse ghost 
resonance may change for surfaces made of different facets (cuts) of the crystal. For instance, at a particular energy a sharp edge 
state on one facet may correspond to a mere ghost resonance when we consider another surface. 

The evolution of the surface-state spectrum in the hybridized system, as described above, makes it clear that that the helical 
characteristic of the surface states will remain unchanged. Strikingly, in some cases the number of surface state branches at each 
edge can change, but the parity of the number of surface-bands will remain odd (see supplementary material Section lA2l where 
the 1 surface state branch gets multiplied to 3). 

Effects of disorder were only briefly and crudely considered here. Nevertheless we can already indicate intriguing features 
which may arise. Following the disordered Fano model for the case of a metallic band overlapping in energy but not in mo- 
mentum, the hybridization will broaden a surface state (say at momentum q) into a ghost, but will also produce sharp (surface) 
energy-eigenstates at energies above and below the energy overlap range, so long as there is an empty patch as a function of 
energy that can support them. An ARPES measurement should be able to observe all these features for metals with an appro- 
priate band structure, where a metallic parasitic band overlaps in one range of energies with topological edge states, but not in 
momentum, and above or below this energy range bulk gaps exist. Tunneling measurements should also be able to observe the 
spectral features we describe here. We will explore this situation more closely in future work. A more mundane prediction of 
the disorder analysis which should be easily accessible in experiment is the quasiparticle lifetime decrease due to increasing 
disorder. 

Topological insulators (and topological phases in general) are characterized by a topological invariant that assumes only 
quantized values, and is thus robust to infinitesimal deformations to the model. The presence of edge states is deeply connected 
to topological order, and the presence of helical edge states in the metallic models we present here would suggest some sort of 
topological order may exist. However, conventional understanding of topological phases depends crucially on having a robust 
global gap in the system. The 'helical metal' phase we consider here is gapless, though in a clean system local gaps in the 
BZ may appear. It is therefore unclear whether one can define a topological invariant in general. However, one incarnation of 
the topological invariant for the topological insulators (in both d = 2, 3), involves Bloch states only at time-reversal invariant 
momenta in the BZ**. If there is a local gap at these points in the BZ, the same topological invariants are still well defined. In the 
'helical metal' phases we introduce here, each of these points can either be locally gapped, or not. It is also unclear what could 
be a topological invariant for the 'helical metal' in the presence of disorder. Going forward, the question whether a topological 
invariant exists for the systems we introduce here is perhaps the most enticing, and we leave its determination to future work. 

The abundance of materials exhibiting topological properties suggests that there must be many materials which are 'helical 
metals' - metals with an odd number of chiral surface states at some energy ranges for each facet. Such materials can presumably 
be found in the vicinity of topological insulators and vice versa, which is supported by the materials observed so far, Sb, 
Bii_3;Sba;, and Bi2 Se^. The spectral effects which we explore here should be accessible in all of these materials and provides 
another challenging system where interaction and disorder may have important and interesting effects. 

It is a pleasure to acknowledge useful discussions with P. A. Lee, O. Motrunich, Z. Hasan, and D. Hsieh. This work was 
supported by the Sherman Fairchild Foundation, by the Packard foundation, Sloan fellowship, and Cottrell fellowship, and by 
the Institute for Quantum Information under NSF grants PHY-0456720 and PHY-0803371. 



Appendix A: Supplementary material 



1. Alternative construction of a helical metal model 



Here we explore an example of an alternative construction of a helical metal model, without introducing new degrees of 
freedom. Both constructions were briefly outlined in Section|IIl We start with the Kane-Mele model'-^, and add second neighbor 
hopping (which involves hopping only on the same sublattice). This realizes the construction leading to a momentum-dependent 
chemical potential. We adjust the second neighbor hopping and overall chemical potential to be strong enough to close the bulk 
gap, while having a near zero value near those momenta at which the surface states appear, on the zigzag edge surface of the 
honeycomb lattice. We explore both this edge as well as the armchair edge. Our model is most simply and succinctly written as 

H = -tiJ2 [4abJc. + ^>'.c]-^lJ24c^''3o. + {a-^b)+ [aL«j73(-^~2^a/3 + «A4/3i/,j-)-h/i.c. + (a^6)] , (Al) 

{ij)a ja {{ij))al3 

where as in Eq. a, b denote the fermion operators on the two sublattices, i,j denote the lattice sites, the spin indices are 
denoted by a, j3, and i^ij is as defined in Refs. 13 (and shown in Fig.|2]l. Repeating the finite strip numerical diagonalization for 
both zigzag and armchair edges, we find the results of Fig. |7(a)| and Fig. |7(b)| demonstrating yet again the presence of helical 
surface states coexisting with a bulk Fermi surface. The surface states in the armchair edge geometry seem as if they have been 
pushed away from overlapping with the bulk states, reminiscent of the exiling effects we discussed in the main text. 
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(a)Zigzag edge 
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FIG. 7: All pictorial conventions are exactly the same as in Fig. |3] including the system size, and the criteria for identifying surface states. 
Here we use parameters A = 0.3, fj. = 0.8, ti = 1 and t2 = 0.4. We show the armchair edge spectrum in |(b)[ and the zigzag spectrum in |(a)| 
The surface states on the two sides of the system perfectly overlap in|(b)| 
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(a)Zigzag edge 
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(b)Armchair edge 



FIG. 8: All pictorial conventions are exactly the same as in Fig. |3] including the system size, and the criteria for identifying surface states. 
Here we use parameters A = 0.08, t2 = —l,t\ = l,h = 2 and = 0.5. We show the spectrum of both zigzag and armchair edge strip 
geometries. The same expulsion of the surface states occurs, as in Fig. [3] and in addition we find the number of surface branches has tripled to 
3 Kramer's pairs on each wall of the strip. 

Note that Some examples of helical metals may also have surface states on some faces but not on others, as would be the 
case in the above model if the second-nearest neighbor hopping were sufficiently strong such that no bulk gap would exist. 
Generically, such surface states will be unstable to strong disorder, and will be subject to a finite life-time as described in Sec. 
llVlabove. 



2. Surface state multiplying 



Another revealing example of the model in Eq. ([T]), shows the helical nature of the surface states being preserved despite 
significant changes in the band structure. The parameters used here in this example are A = 0.08, t2 — —l,ti = l,h = 2 and 
^3 = 0.5, and we, once again, repeat the same numerical calculation on both the zigzag and armchair edged strip geometry. The 
various spectra are displayed in Fig. [8] and show that on both edges, the single surface state branch has multiplied into 3 co- 
propagating surface state branches, albeit not crossing the Fermi energy (though that can be changed with the overall chemical 
potential) so that the total number of helical surface states on one edge has changed from one Kramer's pair of bands to three. 
There are still an odd number of surface state branch pairs, which means the system has the same helical classification. This 
numerical result again confirms the robustness of the helical nature of the surface states in this model, despite strong mixing 
with an ordinary metallic band. We will leave the exploration of the edge-state multiplication effect to future work. 
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